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HIGH  ENERGY  SCALING  C.EN 1 R A L\'/ AT  IONS 
FOR  CO  FLECIRIC  DISCHARGE  LASERS  ( U> 


W.  B.  Lacina  and  G.  L.  McAllister 
Northrop  Research  and  Technology  Center 
Hawthorne,  Califo  nia 

ABSTRACT 

(Unclassified) 

High  optical  extraction  efficiencies  for  CO  electric 
discharge  lasers  have  been  theoretically  predicted  and 
experimentally  demonstrated.  Considerable  interest  for 
military  applications  has  resalted  in  programs  for  develop- 
ment of  both  pulsed  and  cw  CO  EDL's  with  high  power  and 
high  energy  objectives.  In  order  to  make  realistic  perform- 
ance predictions  for  potential  device  design,  it  is  necessary 
to  understand  how  these  devices  scale  in  volume,  pressure, 
and  other  parameters.  In  this  paper,  several  generalizations 
that  result  from  parametric  scaling  of  the  CO  laser  kinetics 
will  be  report  ed  and  compared  with  experimental  data  from 
an  e* -beam  stabilized  device.  Scaling  to  high  energy  and 
high  power  systems  will  be  discussed.  The  limitations  on 
high  pressure  or  line-selected  operation  that  may  be  imposed 
by  resonant  self-absorption  in  the  CO  medium  will  be  dis- 
cussed. Preliminary  investigations  suggest  that  this  effect 
may  be  an  important  pressure-dependent  mechanism,  with 
possible  consequences  for  attainment  of  line  selection 
objectives.  It  will  also  be  shown  that  certain  presently 
unexplained  discrepancies  (relating  to  transient  time  scale, 
spectral  anomalies,  efficiency  degradation,  and  sensitive 
temperature  dependence)  might  all  be  consistently  resolved 
by  this  mechanism. 


1.  INTRODUCTION 

• 

(U)  Because  of  its  favorable  operating  characteristics  of  high  quantum  efficiency  and  high 
specific  output  power,  the  CO  electric  discharge  User  has  attracted  considerate  interest  as  a 
candidate  for  both  pulsed  and  supersonic  cw  devices.  Electron  beam-stabilized  CO  EDL's  have 
been  shown  to  operate'  at  efficiencies  with  specific  energies  >750J// /atr >,  and  high 

output  energies  (>1200  J)  have  been  demonstrated2  with  a 20  liter  pulsed  device.  Five  hundred 
joules  (S00  J)  with  40'ro  efficiency  has  recei.’ly  been  obtained  at  our  laboratory  from  a device 
with  7 liters  ac-  ve  volume,  using  a CO/N2/Ar  = 1 /9/7  mixture  of  190  Torr  at  80°K.  Develop- 
ment of  supersonic  flow  cw  devices  is  currently  in  progress  at  several  laboratories,  and  results 
that  have  been  obtained  at  Mathematical  Sciences  Northwest  and  Northrop  will  be  reported  in  a 
separate  presentation  of  this  Conference. 
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<U)  Development  of  these  devices  for  high  power,  high  energy  military  applications  requires 
that  a number  of  coupled  technological  and  fundamental  physical  questions  be  addressed.  The 
optimistic  theoretical  predictions  are  based  on  laser  kinetic  models  that  make  several  important 
assumptions,  many  of  which  are  necessarily  approximations,  while  others  may  be  omissions.  In 
general,  results  of  calculations  have  agreed  reasonably  well  with  experimental  comparisons, 
although  there  remain  some  discrepancies  which  need  to  be  definitively  resolved  in  order  that 
realistic  predictions  and  device  design  can  be  made  with  confidence.  In  order  to  achieve 
military  objectives  of  h.gh  energy  pulsed  or  „t0„  power  Cw  devices,  the  laws  for  scaling  in 
volume  and  pressure  must  be  correctly  understood.  In  addition,  for  some  applications  for  which 
line  selection  of  only  those  transitions  with  high  atmospheric  transmission  is  required,  the 
consequences  of  scaling  laser  parameters  subject  to  spectral  constraints  must  also  be  assessed. 
Scaling  in  volume  for  these  devices  presents  problems  primarily  on!>  from  component  limitations, 
while  scaling  in  pumping  and  pressure  correspond  to  the  basic  physical  mechanisms  of  the  CO 
laser,  which  are  described  by  the  molecular  and  electron  kinetics  processes  and  the  interaction 
of  the  gain  medium  and  the  optical  extraction. 


(U)  In  this  paper,  several  scaling  generalizations  for  a (pulsed  or  cw)  high  pressure  CO  EDL 
are  presented  and  compared  with  experimental  data  from  a pulsed.  e‘-l,eam  stabilized  device. 

The  scaling  of  parameters  necessary  for  analysis  of  specific  problems  in  terms  of  more  general 
results  >s  discussed,  and  it  is  shown  that  the  transient  evolution  of  the  system  is  mainly  dependent 
only  upon  the  total  energy  deposition  into  CO  as  a function  of  time.  However,  the  validity  of  these 
parametric  scaling  laws  is  based  on  a model  which  may  be  inadequate  to  correctly  describe 
effects  of  high  pressure  (£100  Torr)  operation,  because  of  the  possibility  that  resonant  self- 
absorption may  be  an  important  effect.  All  contemporary  “heoretical  models  of  CO  laser  kinetics 
have  overlooked  the  fact  that  laser  osciUal  on  on  one  line  could  be  absorbed  by  near-resonant 
transitions  of  other  bands.  This  possibility  is  currently  under  investigation,  and  completely  self- 
consistent  calculations  accounting  for  this  phenomenon  will  be  completed  shortly  and  reported  in 
the  near  future.  After  presenting  a summary  of  the  results  of  scaling  the  laser  kinetics  without 
inclusion  of  these  effects,  a discussion  of  several  important  consequences  that  resonant  self- 
absorption may  present  for  high  pressure  devices  will  follow.  It  will  also  be  shown  that  certain 
presently  unexplained  discrepancies  (relating  to  transient  time  scale,  spectral  anomalies, 
efficiency  degradation,  and  temperature  dependence)  might  aU  be  consistently  resolved  by  this 
mechanism.  It  will  be  necessary  to  understand  the  possible  implications  of  this  effect  and  the 
limitations  that  it  may  impose  in  order  to  determine  and  optimize  design  considerations  for  high 
pressure  and  line-selected  devices. 


C.  theoretical  model 

(U)  A molecular  kinetic  model  for  an  elect  ricaUy  excited  (CO,  A r.  He.  . . . ) gas  mixture 
has  been  constructed  for  the  calculation  of  the  radiative  characteristics  of  a CO  laser  oscillator 
or  amplifier.  A computer  program  for  the  numerical  solution  of  the  steady  state  or  transient 
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master  equation  for  the  diatomic  species  has  been  developed  to  predict  vibrational  population 

distributions,  gam  and  satu  ration  pa  rameters , energy  transfer  and  extraction  rates,  conversion 

efficiencies,  output  intensities,  and  spectral  distributions,  The  radiative  calculations  are  self- 

consistent  with  the  saturated  gain  on  all  oscillating  transitions  equal  to  the  losses.  The  plasma 

characteristics,  described  by  a Boltzmann  distribution,  are  adjusted  s elf-consistently  in  the 

calculation  as  a function  of  time  to  account  for  electron  heating  for  a given  temporal  input  power. 

A more  complete  discussion  of  the  assumptions  upon  which  this  model  is  based  has  been  given 

elsewhere,4,5  where  numerical  results  for  a variety  of  operating  conditions  were  presented  and 

compared  with  experimental  data  from  a pulsed  high  pressure  e -beam  stabilized  electric 

discharge  CO  laser  oscillator.  Similar  analytical  work  on  the  CO  EDL  has  been  presented  by 
6-10 

numerous  others 

(U)  It  is  possible  to  make  several  reasonable  approximations  that  allow  certain  general 
scaling  laws  for  a CO  EDL  to  be  formulated.  **  These  generalizations  depend  only  upon  a few 
basic  parameters  which  make  it  possible  to  predict  the  temporal  output  characteristics  of  a 
large  class  of  both  pulsed  and  flowing  cw  CO  systems.  Although  the  present  theoretical  results 
are  limited  to  gas  mixtures  containing  CO  and  monatomic  diluents  only,  comparisons  with 
experimental  data  for  both  CO/Ar  and  CO/N,  are  included.  It  is  found  experimentally  that  the 
excitation  times  required  to  attain  threshold  and  steady  state  agree  with  the  scaling  results 
predicted  for  CO/Ar.  Uncertainties  in  the  CO-N,  and  N^-N^  vv  ra*es  make  theoretical  analysis 
for  CO/N2  mixtures  less  reliable.  However,  the  experimental  results  indicate  that  CO/N^ 
systems  can  be  scaled  similarly. 

(U)  Thus,  good  qualitative  interpretation  of  both  CO/Ar  and  CO /N  data  can  be  made.  It 
can  be  concluded  from  these  results  tha»  a general  CO/X  system  containing  either  monatomic  or 
diluents  can  be  characterized  by  the  temporal  energy  deposition  into  the  CO  vibrational  states. 
Thus,  except  for  the  spectral  details,  it  is  possible  to  design  general  CO/X  systems  to  a very 
good  approximation  without  resorting  to  extensive  computer  calculations. 

(U)  The  physics  of  the  present  theoretical  model  is  described  by  four  basic  systems  of 
equations  that  couple  the  CO  vibrational  state  population  densities,  the  radiation  intensities,  the 
gas  kinetic  temperature,  and  the  plasma  kinetic  parameters  (electron  density  and  temperature) 
as  a function  of  time.  Since  the  general  formulation  of  these  equations  and  the  assumptions  made 
have  been  described  previously,  5 they  will  be  summarized  only  briefly  here.  Certain  approxi- 
mations (most  of  which  are  not  required  for  the  complete  computer  analysis)  make  it  possible  to 
obtain  a variety  of  general  results  in  terms  of  a very  small  set  of  scale  parameters.  In  fact, 
one  basic  parameter  lelated  to  pumping  is  all  that  is  required  to  present  a reasonably  simple 
description  of  the  CO  EDL  under  the  set  of  approximations  to  be  discussed  below. 

(U)  The  physics  of  the  molecular  kinetic  processes  is  described  by  the  master  equation  for 
vibrational  population  densities  n^  of  the  CO  molecule. 
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dn  /dt  * 
v 


< *? R 


VT  (CO.X)  + RVVT  (CO.  CO) 

V V 


V VT  a non  _atim 

4Rvv  1 (CO,  n2i  + rvp  +Rv 


(1) 


with  pumping  and  ration  proce.se.  which  include:  (1)  electron  impact  excitation  of  the 
vibrational  .tate.,  (2)  vibration  io  translation  (VT)  deactivation  proce..e.,  (3)  near  re.onant 
vibration  to  vibration  (VVT)  .ingle  quantum  exchange  collision..  (4)  .pontaneou.  radiative  decay, 
and  (5)  stimulated  emi..ion  and  ab.orption  proces.e..  Source,  for  rate,  have  been  .ummarized 

previ.ouly.  * 

(U)  Radiation  intensities  are  calculated  by  requiring  that  the  gain  coefficient,  for  all 
oscillating  transitions  be  equal  to  the  loi.  coefficient  at  all  time..  The  rotational  level,  are 
assumed  to  be  in  thermal  equilibrium  with  the  molecular  kinetic  temperature,  and  rotational 
cro..- relaxation  i.  a..umed  to  be  .ufficiently  fa.t  that  only  one  rotational  line  can  o.cillate  in 
any  given  vibrational  band.  Lo.se.  .uch  a.  mirror  ab.orption.  output  coupling,  window  lo..e., 
or  intracavity  line  .election  lo..e.  are  all  a..umed  to  be  uniformly  distributed  throughout  the 
la.er  gain  medium.  The.e  cavity  lo..e.  are  a..umed  to  be  .mall  enough  that  the  approximation 
of  a spatially  uniform  intensity  distribution  over  the  length  of  the  active  medium  is  adequate  for 
describing  the  molecular  and  radiative  characteristic,  of  the  medium.  The  optical  resonator  i. 
assumed  to  be  formed  by  two  plane  parallel  mirror,  (one  totally  reflecting,  the  other  with 
reflectivity  R)  separated  by  a distance  L.  From  the  oscillator  gain  conditions.  Rexp[2(o-P)L]  =1, 
the  gain  coefficient  a for  an  oscillating  line  (for  R - 1)  must  be  equal  to  the  threshold  los. 
y=p  + d . R)/2L.  where  pi.  the  loss  coefficient  for  the  cavity  exclusive  of  the  output  coupling. 

(U)  The  analysis  assume,  that  the  plasma  kinetic,  can  be  treated  with  a Boltzmann  electron 
distribution,  with  a self-consistent  adjustment  of  the  electron  density  and  effective  temperature 
a.  a function  of  time  to  produce  a prescribed  temporal  input  excitation  power  characterized  by  a 
rise  and  fall  time.  Physically,  the  choice  of  that  form  for  the  electrical  excitation  power  is 
appropriate  for  an  e*-beam  sustained  CO  device  with  the  rise  time  related  to  the  approach  to 
steady  state  of  the  electron  density  and  the  fall  time  to  the  decay  of  the  su.tainer  voltage  due  to 
the  discharge  of  a capacitor  bank. 


(U)  Gas  heating  from  VV  and  VT  kinetic  processes,  as  well  as  from  rotational  equilibration 
accompanying  stimulated  emission  result,  in  a temporal  rise  in  the  molecular  temperature.  The 
effects  of  temperature  rise  are  included  in  the  analysis. 

(U)  In  order  to  derive  a useful  set  of  generalized  results  that  can  easily  be  scaled,  it  is 
first  necessary  to  simplify  the  equation,  by  introducing  approximation,  that  will  reduce  the  num- 
ber of  scale  parameter,  required  to  a reasonable  minimum  value.  Although  the  additional 
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approximations  to  be  discussed  here  arc  neither  required  nor  invoked  in  the  complete  computer 
analysis,  they  will  be  useful  for  summarizing  the  results  of  the  computer  calculations  in  a 
coherent  way.  The  objective  here  is  to  show  that  a rather  general  class  of  pulsed  CO  oscillators 
can  be  described  quite  accurately  in  terms  of  only  two  basic  scale  parameters  that  define  the 
electrical  pumping  and  cavity  threshold. 

(U)  f or  a high  pressure  oscillator  operating  at  temperatures  typically  in  the  range  60-300°K. 
the  spontaneous  radiation  and  VT  decay  terms  ,n  the  master  equation  can  be  neglected.  Con- 
sistent with  this  approximation,  the  VT  contribution  to  the  kinetic  heating  can  also  be  neglected. 

It  will  be  convenient  to  refer  to  the  initial  values  of  the  temperature  T and  pressures  p 

mol  K KCO’ 

px  in  the  results  of  calculations  to  be  presented  here.  Theoretical  prediction  that  only  a small 
percentage  of  electrical  power  is  converted  to  kinetic  heating  (typically  -10-20%),  even  after 
steady  state  has  been  attained,  have  been  experimentally  verified.  12  Thus,  for  -easonably  dilute 
CC  gas  mixtures,  it  is  a good  approximation  to  neglect  the  temperature  rise  for  pulse  times  that 
are  comparable  to  the  time  required  to  reach  steady  state.  Although  the  calculations  correspond 
to  fixed  number  densities,  the  pressures  will  change  by  only  a small  fraction  during  that  time 
scale.  All  of  the  calculations  to  be  presented  here  correspond  to  case,  for  which  the  electrical 
excitation  was  kept  constant  as  a function  of  time  during  the  pulse,  with  temporal  evoltuion  of  the 

electron  temperature  T . 

c 

(U)  The  coupled  se,  of  equations  of  reference  (5)  with  the  approximation  discussed  above  can 
be  written  in  dimensionless  form  by  introducing  the  following  quantitie.: 

x * n /N 

v v CO 

x = n /N 

e e CO 

T = pcoT 

' * 1 4 


i = A v /A  v 

CO  tot 


pCO<  Vrel(C°’  C°l>-°Pt(CO.CO) 
2 Px<  vrel<CO,  X)  >T0pt(CO,  X) 


\ * VPCO 

f = I{/pco 

E * E/p ™ 


where  xy  are  the  relative  populations,  xg  is  a measure  of  the  degree  of  ionization,  f is  related  to 
the  number  of  thermal  degrees  of  freedom,  and  % is  the  fractional  percentage  of  total  optical 
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line  broadening  caused  by  CO.  Y,  W,  I,  and  E are  scaled  threshold  loss  (cm*1),  power  density 
(w/cm  /Torr  (CO)],  intensity  [ W/cm2 /To rr2(CO|],  and  energy  density  [j/cm3/Torr(CO)].  The 
definition  of  the  parameter  4 given  here  applies  to  a situation  in  whi*  pressures  are  high  enough 
(220  Torr)  that  optical  lines  are  predominantly  pressure  broadened. 

(U)  If  these  scale  parameter,  are  introduced,  the  master  equation  expresses  the  evolution 
of  the  relative  populates  xy  as  a function  of  scale  time  r = p^t.  This  U to  be  expected  since 
the  electrical  excitation  is  kept  constant,  and  the  dominant  kinetic  processes  are  (two  body)  VV 
collisions.  Thus,  the  physics  of  the  molecular  vibration  states  depends  explicitly  only  upon  the 
kinetic  temperature  and  the  scaled  electrical  excitation  W#  = V^/p^.  Implicitly,  however,  it 
depend,  also  upon  the  scaled  threshold  loss  coefficient  Y for  the'  cavity  due  to  the  oscillation 
condition,  and  upon  the  parameter  f due  to  the  kinetic  heating. 

<U)  I,  is  apparent  that  the  s-t  of  parameter.  (T^,  Y . f)  completely  characterize  the 
CO  EDL  under  these  approximations;  ah  other  input,  to  the  mathematical  model  consist  only  of 
fundamental  physical  constants.  In  the  following  section,  a variety  of  numerical  result,  based  on 
the  complete  computer  analysis  for  this  model  will  be  presented  and  discussed  to  show  the 
relative  sensitivity  of  these  parameters. 


3.  THEORETICAL  RESULTS 

(U)  Computer  calculation.  Here  carried  out  for  a fixed  value  f * 16  (corresponding  to  a 
dilute.  10T.  CO  mixture).  The  kinetic  temperatures  T^  included  (60.  100,  150,  200.  300)°K, 
the  scaled  electrical  pumping  rates  (0.2,  0.5,  1.0,  2.0,  5.0,  10.  20.  50,  100,  200)W/Cm3/TorrZ, 
and  the  scaled  threshold  loss  coefficients  (0.1,  0.2,  0.5,  1.0,  3.0,  7.0,  10,  20)%/cm.  From 
the  result,  of  any  given  case  specified  completely  by  the  set  of  four  parameters  (T  , W , Y , f), 
the  general  class  of  problems  that  can  be  scaled  exactly  from  those  result,  is  quit  eTestrictive; 
it  allows  only ^he  variation  of  the  CO  partial  pressure,  provided  that  the  electrical  excitation  W 
is  scaled  by  pco,  and  that  the  threshold  level  and  gas  mixture  ratio  are  held  constant.  However, 
except  for  the  scaled  excitation  power  W , most  of  these  parameter,  do  not  have  a sensitive 
effect  on  the  power  response  of  the  system  as  a function  of  time  (although  detailed  spectral 
prediction,  may  have  a more  critical  dependence).  Thus,  the  class  of  problem,  that  can  be 
analyzed  from  a given  case  is  often  more  extensive  than  that  jus,  described.  The  reasons  for 
this  are  the  following  (1)  the  most  important  parameter,  the  scaled  excitation  is  a measure 
of  the  electrical  pumping  per  molecule  relative  to  the  rate  of  VV  cross- relaxation'  and  thus 
specifies  the  ratio  of  the  two  dominant  kinetic  processes  in  the  CO  EDL.  (2)  A pulsed  high 
pressure  CO  EDL  typically  operates  in  a regime  where  oscillation  is  high  above  threshold,  and 
thus,  results  for  optical  power  extraction  will  usually  be  insensitive  to  the  scaled  cavity  loss 
coefficient  Y.  Unless  the  ratio  'H J Y becomes  very  small  the  medium  intensity  will  be  high 
above  the  saturation  intensity,  and  the  optical  output  will  be  independent  of  Y.  (3)  The  tempera- 
ture rise  for  pulse  times  comparable  to  the  time  required  to  , each  steady  state  is  small,  so  the 
result,  do  not  depend  critically  upon  f,  providing  it  is  large.  This  parameter  is  mainly 
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important  for  detailed  predictions  of  spectral  distribution.  (4)  Results  of  calculations  for 
systems  operating  with  parameters  in  the  above  ranges  do  not  show  a significant  sensitivity  to 
depend enc e on  T 

mol 


(U)  Figure  1 shows  plots  of  the  quantum  power  efficiency  (for  the  total  optical  conversion  of 
electrical  input)  for  a variety  of  different  values  of  W as  a function  of  p t,  for  T = 100°K 
Results  of  calculations  show  that  the  most  important  parameter  for  small  values  of  Y is  the 
scaled  excitation  power  W^,  and  that  the  optical  output  is  relatively  insensitive  to  Y for  small 
values  of  Y . It  is  apparent  from  Figure  1 (and  from  similar  results  for  other  values  of  Y and 
Tmol)  that  the  optical  outPut  efficiency  displays  a remarkably  similar  behavior  over  a wide 
range  of  the  parameter  W^.  In  fact,  these  results  can  be  used  to  show  that,  for  a given  tempera- 
ture  Tmol  and  a fixed  value  of  quantum  efficiency  n/na  (where  n<a  is  the  steady  state  value),  a 
loglog  plot  of  pCQt  versus  the  scaled  power  is  approximately  linear.  For  example,  Figure  2 
shows  such  a plot  for  the  turn-on  time  PCC)to  (obtained  by  extrapolating  the  curves  of  Figure  1 to 
the  origin  and  neglecting  a small  initial  bump  that  often  occurs  in  the  calculations)  versus  W 

0 

The  slope  of  these  curves  is  approximately  -5/6  and  is  independent  of  T and  V;  the  variation 

mol 

in  the  magnitude  of  the  curves  is  less  than  a factor  of  two  within  the  range  of  temperatures 

chosen.  The  turn-on  time  corresponds  to  the  point  where  n = 0;  however,  it  is  apparent  that  the 

curves  of  Figure  1 remain  similar  as  time  evolves.  Thus,  for  a fixed  temperature  T the 

mol' 

timet  required  to  attain  any  specified  value  h/q  for  the  optical  extraction  efficiency  satisfies 


,?,5/6 

e = Pco1  e * constant, 


(3) 


which  implies  that,  within  a range  of  values  for  Y which  are  suitably  small,  a "universal"  plot  of 
efficiency  versus  the  parameter  e = p^tW^^  can  be  constructed.  As  an  example,  all  of  the 
plots  for  the  conditions  of  Figure  1 can  be  summarized  as  shown  in  Figure  3,  (The  value  of  the 
steady  state  efficiency  n,,  depends  only  slightly  upon  and  V . ) The  significance  of  the 

parameter  e can  be  seen  by  noting  that  if  the  exponent  of  W in  Eq.  (3)  were  unity  instead  of 
5/6,  e would  represent  the  total  electrical  energy  E = /pCQ[ J // /Torr(CO)]  which  has  been 

deposited  into  CO  vibrational  energy  up  to  time  t.  Therefore,  the  response  of  the  system  is 
approximately  a functior  only  of  the  specific  input  energy  E(which  is,  of  course,  a function  of 
time  and  the  excitation  rate).  The  power  efficiency  can  be  integrated  to  give  a universal  curve 

for  the  energy  efficiency  «/nx  as  a function  of  the  parameter  e,  and  this  result  is  also  presented 
in  Figure  3. 


(U)  Among  the  important  conclusions  to  be  drawn  from  the  above  remarks  are  the  following. 
Characteristic  times,  such  as  turn-on  time  or  time  required  to  reach  steady  state,  satisfy  a 
relation  given  by  Eq.  (3).  For  the  range  of  temperatures  and  for  the  reasonably  small  values  of 
Y discussed  here,  turn-on  occurs  after  a specific  electrical  energy  E - 0.  5 - 1. 0 J/f/Torr  (CO 
has  been  deposited,  and  steady  state  is  attained  after  an  energy  Egg  - 1.7  - 2.6  J///Torr  (CO  has 
been  deposited.  (These  numbers  correspond  to  a specific  pumping  rate  W = 1,0  W/cm3/Torr 
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2 - —1/6 
(CO)  , for  other  values  of  W , Eq.  (3)  shows  that  these  deposition  energies  increase  as  W ). 

C * e -1/6 
The  typical  range  of  values  for  W in  the  data  to  be  presented  later  is  sufficiently  small  that 

variation  is  not  significant,  and  thus,  experimental  results  will  be  plotted  in  terms  of  the  more 

physically  significant  quantity  E = W^t/p , 


(U)  It  is  interesting  to  characterize  the  CO  medium  by  a saturation  intensity.  Although 
the  full  kinetic  model  contains  simultaneous  oscillation  on  several  lines,  it  is  convenient  to 
neglect  this  feature  for  the  moment,  and  to  investigate  how  the  gain  saturates  as  a function  of  the 
total  intensity  in  the  medium.  F rom  the  analysis  of  a simple  two-level  model,  it  can  be  shown 
that  the  saturation  gain  a (I)  satisfies 


o(I)  = o(l+  I/I  f 1 , (4) 

o s 

where  the  saturation  intensity  1^  is  independent  of  the  electrical  pumping.  We  shall  show  that, 
at  steady  state,  the  CO  EDL  follows  a simple  dependence  of  this  form,  and  an  estimate  of  the 
scaled  saturation  intensity  as  a function  of  T wiH  be  derived.  In  terms  of  scaled  quantities, 
Eq.  (4)  becomes,  upon  setting  o(I)  = Y , 

Y + q w /f  = o , (5) 

® e s o 

where  the  steady  state  cc  ;ion  Yf  = q flr  has  been  used.  Thus,  for  a fixed  value  W , the 
’ ® e e 

steady  state  oscillator  power  efficiency  q^and  the  cavity  threshold  Y are  linearly  related  if  the 
laser  gain  medium  saturates  according  to  the  relation  <-f  Eq.  (4).  That  is,  if  q (0)  is  the  limiting 
value  of  the  steady  state  efficiency  for  Y—0  (lossless  cavity),  then 


Y = -lq  -n<0)]W  /l  (6) 

® ® e s 

(U)  Figure  4 shows  q plotted  versus  Y for  T = 100°K,  for  two  different  values  of  W , 

® _ mol  _ e 

and  as  Eq.  (6)  shows,  the  slope  | dY/dq^l  of  these  curves  is  proportional  to  W Furthermore, 
the  proportionalitv  constant  determines  the  scaled  saturation  intensity,  which  is  determined  from 
Figure  4 to  be  I * 0.8  W /cm  /Torr  . Similar  plots  for  a fixed  value  of  W and  variable  tempera- 
ture  are  shown  in  Figure  5.  All  of  these  curves  are  linear  to  a very  good  approximation, 

and  the  saturation  intensities  obtained  from  their  slopes  are  summarized  in  Table  I.  During  the 
transient  evolution,  the  relation  between  the  quantum  efficiency  and  the  cavity  threshold  does  not 
seem  to  follow  a linear  relationship. 


TABLE  I.  SCALED  SATURATION  INTENSITY  (U) 
(Unclassified) 


Tm.l'°KI 

60 

100 

150 

200 

300 

2 2 

Is(W/cm  /Torr  ) 

0.6 

0.  8 

1.6 

2.5 

6.  4 
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Small  signal  gain  characteristics  can  be  obtained  from  Eq.  (5)  with  9=0.  For  values  of  pumping 

and  cavity  threshold  considered  here,  it  is  apparent  that  the  small  signal  gain  is  proportional  to 

W , 
e 


S = n (0)  W /I  (7) 

o ® e s 

For  efficient  operation  of  the  oscillator  (high  above  threshold),  the  necessary  condition  is 

a »Y  and  therefore  (since  q -1), 
o ® 

w /Y»r  . (8) 

e s 

4.  EXPERIMENTAL  COMPARISONS  AND  CONCLUSIONS 

(U)  The  electrical  excitation  times  required  to  reach  laser  threshold  and  steady  state  have 
been  measured  for  CO/Ar  and  CO/N^  gas  mixtures  and  are  compared  below  to  the  scaling 
generalizations  described  in  Section  3.  The  laser  used  for  these  measurements  was  an  ("-beam 
stabilized,  electric  discharge  CO  laser  with  a 10  cm  x 10  cm  x 100  cm  discharge  region.  The 
operating  temperature  was  80°K  with  gas  pressures  typically  in  the  range  from  50  - 200  Torr. 

The  nominal  operating  parameters  were  E/N  : 1,1  x 10  cm^,  sustainer  current  density 

2 - 2 
1^  = 2A/cm  , and  e -beam  current  density  1^  = 25  mA/cm  . Output  energies  exceeding  500J/ 

pulse  have  been  obtained.  The  electrical  excitation  pulse  length  was  75  ps  and  the  optical  pulse 

length  varied  between  20  - 100  ps  depending  on  the  electrical  excitation  rate  and  the  gas  mixture 

used. 

(U)  The  laser  output  energy  was  measured  with  a ballistic  thermopile  calorimeter  (Hadron 
Series  117)  and  the  temporal  pulse  shape  was  monitored  with  a Ge:Au  detector.  Figure  6 
illustrates  a typical  pulse  shape  for  CO/N2  with  a nearly  constant  electrical  excitation  rate  for 
75  ps.  In  this  example  laser  threshold  is  reached  38  ps  after  the  electrical  excitation  is  initiated 
and  the  output  lasts  considerably  longer  than  the  excitation  pulse  length.  This  long  decay  is 
typical  of  CO/N^  gas  mixes  due  to  the  time  required  for  energy  stored  in  the  vibrational  levels  of 
nitrogen  to  be  transferred  back  to  the  CO  molecules.  The  pulse  shapes  for  CO/Ar  mixtures  look 
similar  except  that  the  pulse  decay  time  is  shorter,  since  no  energy  is  stored  in  the  argon. 

(U)  The  laser  threshold  time  t is  taken  to  be  that  time  from  initiation  of  the  sustainer 

o 

current  to  the  point  where  the  laser  intensity  of  3 - 5%  of  the  maximum  intensity  as  indicated  by 

the  arrow  in  Figure  6.  This  threshold  is  generally  sharp  and  the  time  is  therefore  well  defined. 

The  steady  state  time  t is  defined  as  that  time  required  to  reach  -95%  of  the  maximum  value, 

also  indicated  by  an  arrow  in  Figure  6.  This  value  is  not  as  sharply  defined  as  t and  consider- 

o 

able  spread  in  the  data  may  occur  due  to  the  interpretation  of  when  steady  state  occurs. 

(U)  It  was  shown  in  Section  3 that  the  inverse  of  the  threshold  and  steady  state  times  for 

CO/Ar  laser  gas  mixtures  should  scale  as  p W5^6,  where  W is  the  electrical  excitation  rate 

vU  c c 
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[W/cm  /Torr  (CO)].  For  the  range  of  variables  considered,  the  exponent  (S/6)  can  be  taken  to  be 
unity  with  the  resulting  physical  interpretation  that,  regardless  of  the  eacitation  rate,  fixed 
specific  energies  [ J/ / /Torr(CO) j must  be  deposited  to  attain  threshold  and  steady  state.  Experi- 
mental values  of  threshold  and  steady  state  times  were  measured  for  CO/Ar  gas  mixtures  from 
76  to  ISO  Torr.  The  inverse  times  are  plotted  in  Figure  7 and,  as  predicted  by  the  scaling 
generalizations,  are  proportional  to  W /P(-q-  From  Figure  7 it  is  concludeu  that  approximately 
0.9  J///Torr  (CO ) must  be  deposit  ed  in  order  to  reach  threshold,  which  is  in  good  agreement 
with  the  scaling  generalizations,  which  predict  that  0.S  - 1.0  J/  / /Torr  (CO  is  necessary.  The 
experimental  value  for  the  deposited  energy  required  to  attain  steady  stat r , -2.0  J///Torr  (CO), 
is  also  in  good  agreement  with  the  scaling  predictions  of  1 . 7 - 2.6  J///Torr  (CO). 

(U)  The  calculations  discussed  in  Section  3 are  valid  only  for  gas  mixtures  of  CO  and 
monatomic  diluents.  However,  the  scaling  generalizations  appear  to  be  valid  for  CO/diatomic 
systems,  provided  that  only  the  fractional  energy  into  the  CO  system  is  used  for  scaling. 
Experimental  measurements  of  the  threshold  and  steady  state  times  for  a fixed  CO/N^  (16) 
mixture  are  shown  in  Figure  8.  It  is  apparent  from  this  graph  that  the  inverses  of  these  times 
scale  wdth  W just  as  they  do  for  CO/Ar  mixtures  although  the  required  total  energy  deposited  is 
different.  This  is  to  be  expected  because  the  deposited  energy  is  split  between  the  CO  and  N^. 

It  would  seem,  in  fact,  that  the  required  energy  deposited  in  the  CO  vibrational  system  should  be 
the  same  regardless  of  the  diluent  and  that  the  fractional  amount  deposited  into  nitrogen  can  be 
determined  from  the  experimental  data.  That  is,  a specific  energy  EfCOl/p^Q  must  be  deposited 
into  vibrational  excitation  of  CO  in  order  to  attain  any  given  efficiency. 

(U)  Additional  evidence  is  provided  for  the  previous  supposition  by  measurements  of  the 

threshold  time  dependence  on  the  fractional  content  of  CO,  f__  = p„_/p  . This  data,  shown  in 

CO  rCO  rtot 

Figure  9,  indicates  that  the  total  specific  energy  required  to  achieve  threshold  increases  linearly 
with  f which,  when  extrapolated  to  f^Q  = 1,  gives  - 0.  7J// /Torr(CO),  whic)\  agrees  approxi- 
mately with  the  previous  results  for  CO/Ar  mixtures.  If  it  is  assumed  that  the  required  threshold 
energy  deposited  in  CO  is  always  0.  7J //  /Torr(CO),  the  remaining  fraction  into  nitrogen  increases 
linearly  with  fv  = 1 - f , as  is  also  shown  in  Figure  9.  That  is, 


E (tot) /p 

o lOt 


= f E 
CO  o 


(COI/Pco 


+ f 


n2 


(8) 


From  the  experimental  value  for  E (CO)/p  , it  is  determined  that  E (N  )/p  = 0.  13J///Torr(N  ). 

o CO  o 2 N 2 

The  physical  interpretation  of  the  experimental  data  presented  in  Figure  9 is  tnat  the  energy 
deposition  per  molecule  is  six  times  higher  than  CO  than  for  N^.  Thus,  it  appears  that  the 
scaling  generalizations  derived  for  CO/Ar  mixtures  can  also  be  applied  to  systems  with  diatomic 
diluents,  and  that  accurate  knowledge  of  the  large  number  of  kinetic  rates  required  for  a complete 
analysis  of  the  CO/N^  system  is  not  essential  for  prediction  of  approximate  laser  performance. 
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5.  RESONANT  SELF  ABSORPTION 

(U)  An  important  pressure  effect  that  ii  not  included  in  the  physics  of  the  present  kinetic 
model  described  in  Section  2 is  the  self-absorption  (or  enhanced  gain)  of  laser  radiation  in  the 
medium  by  coincidental  resonances  originating  from  other  P and  (mostly)  R transitions.  As  the 
pressure  and  temperature  increase,  the  collision  broadened  half-widths  of  all  the  CO  lines 
increase,  resulting  in  several  significant  overlaps  for  almost  every  transition.  Thus,  it  is 
likely  that  the  simple  pressure  dependence  described  ir.  the  previous  sections  for  scaling  tlwCC 
laser  is  probably  not  completely  accurate,  and  that  such  a model  may  be  seriously  invalid  for 
prediction  of  detailed  spectral  distributions.  The  impact  of  this  phenomenon  on  previous  com- 
parisons between  theory  and  experiment  will  be  described  below,  and  a variety  of  discrepancies 
that  are  currently  unexplained  will  be  shown  to  have  a consistent  resolution  based  on  the  hypothesis 
of  resonant  self-absorption.  Initial  calculations  have  been  carried  out  using  the  existing  kinetics 
code  for  an  atmospheric  pressure  CO  EDL  at  1C0°K,  pu  nped  at  rates  typical  of  present  experi- 
mental devices  (10  kW/cm  ),  and  these  show  that  self-absorption  competes  strongly  with  all 
other  kinetic  and  excitation  rates.  A refinement  and  extension  of  the  laser  kinetics  computer 
code  has  been  undertaken  to  properly  include,  in  a completely  self-consistent  way,  all  of  the  gain 
and  absorption  processes  to  determine  radiation  intensities  for  an  oscillator.  These  results  will 
be  available  shortly,  and  will  be  presented  in  the  near  future. 

(U)  The  importait  consequences  that  are  anticipated  from  this  analysis,  assuming  that  the 
effects  are  indeed  significant,  will  be  described  briefly  here. 

5.  1 TRANSIENT  TIME  SCALE 

(U)  The  absorption  of  the  medium  is  expected  to  increase  laser  threshold  time,  and 
probably  to  lengthen  the  transient  evolution  to  steady  state.  Experimental  comparisons  made 
previously,  5 as  well  as  those  reported  above  in  Section  4,  are  consistent  with  this  possibility, 
since  experimental  times  for  onset  of  laser  oscillation  were  always  greater  than  theoretically 
predicted. 

5.2  SPECTRAL  ANOMALIES 

(U)  Theoretical  calculations  consistently  show  that  predicted  output  spectral  distributions 
should  lie  about  one  or  two  v-hands  lower  than  what  is  usually  observed.  Furthermoi  -,  experi- 
mental output  spectra  often  show  that  certain  transitions  are  always  missing,  and  it  was  from  this 
observation  that  the  hypothesis  of  resonant  self-absorption  by  accidental  coincidences  among  the 
CO  transitions  was  first  conceived  to  be  an  important  mechanism.  Typically,  the  resonances 
result  in  one  of  two  possible  effects:  1)  The  absorptions  are  so  strong  that  laser  oscillation  is 

permanently  inhibited  from  certain  transitions,  or  2)  oscillation  occurs,  but  output  coupling 
efficiency  is  degraded  and  higher  vibrational  levels  arc  absorptively  pumped.  The  first  effect 
can  have  serioi*  consequences  for  line  selection,  since  it  becomes  possible  that  certain  good 
atmospheric  transmission  lines  le.  g.  , 6—5,  P(10)]  may  never  be  attainable  from  a high  pressure 
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device.  (6  — S,  P(10)  has  never  been  observed  from  oar  e -beam  device  even  at  total  pressures 
-85  Torr  at  100°K.)  The  second  possibility  may  have  serious  con«equences  for  degradation  of 
the  optical  extraction  efficiency,  gas  heating,  and  temperature  dependence  of  the  predicted 
characteristics.  These  are  discussed  in  the  next  section.  Table  II  presents,  as  an  example,  a 
summary  of  lines  which  are  resonant  with  transitions  which  are  desirable  for  atmospheric 


t ransmis  sion. 


TABLE!  II.  RESONANT  TRANSITIONS  FOR  SELECTED  CO  LINES  WITH 
GOOD  ATMOSPHERIC  TRANSMISSION  (Ul 
(Unclas  sified) 


Laser  Line 


E requenc y 
(cm**  I 


1 977. 264 
1 973 . 285 


Attenuation 
Range  (km! 


Resonant 

Transitions 


-1 

^v(cm  ) 


4 - 3 P(9) 


2003. 154 


1982. 754 


1978. 575 
1974. 362 
1957. 189 


2037. 113 
2033. 132 
2029.  1 17 


7- 6  R(4) 

8- 7  R(  1 2 ) 

4-3  P(20) 

9- 8  R(  1 4 ) 


9-8  R(6 ) 


6-5 


6- 5 

7- 6 


4 — 3 P(1Q) 


2025. 068 


2012. 723 


4-3  P(  1 5 ) 2004. 326 

4-3  P(20 ) 1982. 765 


3-2  P(  5)  2071.  146 


3 - 2 P(8)  2059.  200 


3-2  P(10)  2051. 066 

3 - 2 P(13)  2038.615 


3-2  P(15)  I 2030. 148 


6-5 


9 — 8 RIM) 

4-3  R(l) 
6-5  R(  1 7) 


2 - 1 P(  1 9) 

7—6  R(  15) 


Calculations  of  D.  K.  Rice  based  on  molecular  absorption  through  horizontal  path 
at  sea  level  in  a mid  latitude  Winter  atmosphere.  This  model  utilizes  the  McClatchey 
AFCPi^  major  atmospheric  constituent  data  covering  wavelength  from  1 pm  to  20pm. 
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5.  3 OPTICAL  EXTRACTION  AND  EFFICIENCY 

(U)  Theoretical  calculations  of  bulk,  optical  extraction  have  generally  shown  that  quantum 
efficiencies  seem  to  approach  a universal  steady  state  value  of  approximately  90%,  with  littie 
dependence  on  molecular  temperature  or  electrical  pumping  rates.  For  example,  the  family  of 
curves  in  Figure  1 show  littie  sensitivity  to  the  (scaled)  pump  rate  W /P^  , and  similar  curves 
are  predicted  at  high  temperature,  suggesting  that  there  is  a negligible  dependence  on  tempera- 
ture. It  is  found  experimentally,  however,  that  room  temperature  operation  is  significantly 
less  optimistic  than  that  predicted  by  existing  theoretical  models,  and  quantum  efficiencies  of 
-90%  have  not  yet  been  attained  even  at  low  temperature.  Degradation  at  optical  extraction 
efficiency,  as  well  as  the  poor  agreement  with  high  temperature  predictions,  can  both  be  quali- 
tatively explained  by  the  resonant  absorption  hypothesis.  The  optical  extraction  per  unit  volume 
is  determined  by  a combination  of  quantum  efficiency  and  output  coupling  efficiency,  where  the 
latter  is  the  ratio  of  the  output  coupling  loss/pass,  divided  by  the  totai  optical  loss/pass.  If  a 
line  has  attained  sufficient  gain  to  lase,  absorption  in  the  medium  is  effectively  added  to  other 
losses  (e.  g.  windows)  in  the  cavity,  and  thus,  absorption  in  the  medium  (which  is  typically  in 
higher  v-bands)  will  resuit  in  a degraded  coupling  ratio.  Furthermore,  the  absorptive  pumping 
of  higher  vibrationai  levels  wiil  degrade  the  predicted  quantum  efficiency,  since  VV  heating  wiil 
be  increased  and  the  capacity  to  extract  energy  in  the  form  of  laser  radiation  is  reduced. 

5.4  TEMPERATURE  DEPENDENCE 

(U)  At  higher  temperatures,  tor  a given  density  the  collision  broadened  half-width  is 
increased,  and  therefore  the  number  of  resonant  absorption  processes  that  may  contribute  is 
increased.  Furthermore,  the  P(J)  transitions  are  shifted  to  higher  J values,  which  may  result  in 
R(J)  resonances  with  more  severe  absorptive  characteristics.  Lastly,  gain  and  absorption  are 
strong  functions  of  temperature,  and  although  the  oscillator  operatec  with  saturated  gain  equal  to 
loss,  the  increased  absorptive  pumping  of  higher  vibrational  ievels  may  explain  the  poor  perform- 
ance characteristics  of  high  pressure  room  temperature  devices.  It  has  often  been  speculated 
that  the  unexplained  discrepancy  in  temperature  dependence  may  be  due  to  inadequate  knowledge 
of  the  VV  cross-relaxation  rates.  However,  it  can  be  shown  by  a simple  argument  that,  if  the 
functional  dependence  of  the  VV  rates  on  the  ievels  v is  held  fixed,  adjusting  the  magnitude  of  the 

VV  rate  constant  cannot  account  for  the  insensitivity  that  the  quantum  efficiency  shows  to  tempera- 

o 

ture.  To  see  this,  suppose  that  the  family  of  curves  shown  in  Figure  1 corresponded  to  T = 300  K; 

then,  for  fixed  excitation  W , the  effect  of  changing  the  VV  rate  can  be  determined  by  examining 

e 2 

curves  for  different  values  of  W /p  since  that  parameter  is  merely  a measure  of  the  eiectricai 

C LU 

pumping  per  molecule  relative  to  the  rate  of  VV  cross  relaxation.  Since  the  family  of  curves  in 
Figure  1 span  more  than  two  orders  of  magnitude  for  the  parameter  ^/p^Q,  i*  *•  not  likely  that 
better  agreement  between  theory  and  experiment  wiil  be  attainable  by  ad>isting  the  VV  rates  with 
only  a multiplicative  constant.  Other  possibilities,  for  explanation  of  th«  anomolous  temperature 
dependence  are  that  the  structure  of  the  VV  rate  matrix  is  drastically  different,  or  that  VT 
processes  from  high  vibration  ievels  are  much  faster  than  predicted  from  extrapolation  of  SSH 
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theory  (fit  to  1 —0  experimental  data),  which  is  not  likely.  Resonant  self  absorption  may,  therefore, 
resolve  the  unexplained  sensitivity  to  temperature  dependence. 


S.S  FREQUENCY  PULLING  AND  PUSHING 

(U)  Because  the  gains  for  all  of  the  CO  transitions  are  now  sums  of  Lorentzians  for  all  of 

the  resonant  lines,  the  frequency  dependence  of  the  (tain  is  distorted.  It  is.  therefore,  possible 

that  one  of  the  important  consequences  of  s elf -abso rpt ion  may  be  the  shifting  of  the  oscillation 

frequency  from  the  line  center.  Figure  10  illustrates  a case  for  a line  with  gain  centered  at 

frequency  v with  a neighboring  absorbing  line  located  at  ( c • where  iu  is  the  collision 

° o 

broadened  half-width.  The  net  gain,  as  a function  of  frequency,  is  a sum  of  these  two  Lorentzians 

and  is  also  shown  The  line  center  for  the  net  gain  occurs  at  a frequency  shifted  away  from  u 

0 

(in  this  case  pushed)  by  an  amount  4c  . Similar  effects  resulting  in  pulling  can  result  if  the 
neighboring  resonant  line  exhibits  gain  rather  than  absorption.  The  magnitude  of  these  frequency 
shifts  may  be  important  for  atmospheric  transmission.  If  they  are  important  (e.  g.  >.02  cm*') 
the  effect  will  not  be  detected  by  experimental  absorption  measurements  based  on  low  pressure 
CO  laser  probes  oscillating  on  transition  line  centers. 
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THRESHOLD  TIME  pC()t0  (ps- Torr) 


EXCITATION  We/P^0  (W/cm3/Torr2) 


FIGURE  2 SCALED  THRFSHOLD  TIME  pCQt0  AS  A FUNCTION  OF  SCALED  EXCITATION 
WjpcoZ  FOR  SEVERAL  VALUES  OF  LOSS  Y AT  100°K.  (U) 

(Unclassified  figure) 
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QUANTUM  POWER  ANO  ENERGY  CONVERSION  EFFICIENCIES 
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FIGURE  3.  GENERALIZED  POWER  AND  ENERGY  CONVERSION  EFFICIENCIES  AS  A 
FUNCTION  OF  Pcot(  We/PCOZ,5/6  AT  l00°K-  (Nole  th,tl  witl,out  the  5/6  exponent,  the 
unit*  of  the  horizontal  axi*  would  corre»pond  to  energy  depo»ition  in  J/ / / Torr(CO))  (U) 

(Uncial* l fie d figure) 
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FIGURE  4 STEADY  STATE  POWER  EFFICIENCY  AS  A FUNCTICJNOF  SCALED  THRESHOLD 
LOSS  COEFFICIENT  S'  FOR  AN  OSCILLATOR,  FOR  TWO  VALUESIF  SCALED  EXCITATION 
We  s we/PcO^-  Note  that  the  slopes  ol  these  curves  are  proportions!  to  I /W  where  Tf 
is  the  scaled  saturation  intensity.  (U) 


(Unclassified  figure) 
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FIXTURE.  AS  A FUNCTION  OF  We/pc0.  (U) 


(UnclaisLfLed  figure) 
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SPECIFIC  ENERGY  (J///TORR) 
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FIGURE  9.  EXPERIMENTAL  SPECIFIC  ENERGIES  AS  A FUNCTION  OF 
FRACTIONAL  CO  CONTENT  REQUIRED  FOR  ATTAINMENT  OF  LASER 
THRESHOLD  FOR  CO/Nz  GAS  MIXTURES.  (U) 


(Unclaaaified  figure) 
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FIGURE  10.  FREQUENCY  DEPEN  DENT  CAIN  FOR  OVERLAPPING 
COLLISION  ALLY  BROADENED  NEAR  RESONANT  TRANSITIONS.  (U) 

(UncUsaifir.d  figure) 
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